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ABSTRACT 

We have conducted a variability survey of the local group galaxy M33 using g',r', and i' 
observations from 27 nights spanning 17 months made with the MegaPrime/MegaCam 
instrument on the 3.6 m CFHT telescope. We identify more than 36000 variable sources 



with g', 



< 



24 out of approximately 2 million point sources in a one square degree 



field of view. This increases the number of known variables in this galaxy by more 
than a factor of 20. In this paper we provide a brief description of the data and a 
general overview of the variable star population which includes more than 800 candi- 
date variable blue and red supergiant stars, more than 2000 Cepheids, and more than 
19000 long period variable AGB and RGB stars. 
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1 INTRODUCTION 

The study of variable stars in local group galaxies has played 
an important role in our understanding of the universe. Tra- 
ditionally variable stars have served as one of the prime tools 
for determining the local distance scale. The discovery by 
Leavitt (1908; Leavitt and Pickering 1912) of the period- 
luminosity relation for Cepheids in the Magellanic clouds 
and the use of this relation by Hubble (1925, 1926, 1929) on 
Cepheids in NGC 6822, M33 and M31 enabled the modern 
view that these are all galaxies located outside the Milky 
Way. 

In recent years large-scale photometric variability sur- 
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veys have revolutionized the st udy of variable stars . Dark 
matter searches such as EROS (lAubourg et alJll993lh MA- 
CHO dAlcock et ai1ll997Tl. MOA <Noda et all l2002f) and 
OGLE ifudldslri'^^ibiak and SzvmahskiT^97i) have re- 
leased catalogues of more than two hundred thousand vari- 
ables in the Galactic Bulge and more th an sixty thou- 
sand variables in the M agellanic Clouds (Wozni ak et alJ 
12001 IZebrun et alJlioOlTl . These surveys have shown that 
completely new science can be done with very large sam- 
ples of variable stars. For example, a few of the re- 
sults based on these catalogues include the discovery of 
multiple period-lum inosity relations for Mira-like variables 
JWood et alJll999l) the d iscovery of Beat Cepheids in the 
LMC jAlcock et aljll995T) and the discovery of a period- 
luminosity relati on for ellipsoidal va riations observed in 
LMC red giants JSoszvnski et al.ll20o4) . 

While the Magellanic Clouds have been surveyed exten- 
sively by the microlensing searches, larger and more distant 
galaxies such as M31 and M33 have only been targeted with 
CCDs by relatively narrow-field variability surveys. With 
the advent of wide-field mosaic imagers it has become pos- 
sible to perform a complete survey of these galaxies with- 
out the use of a dedicated telescope. Motivated by this op- 
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portunity, we have conducted a variability survey of the 
nearby galaxy M33 (01:33:51.00 +30:39:36.7 J2000) using 
the MegaPrime/MegaCam instrument on the 3.6 m Canada- 
France-Hawaii- Telescope (CFHT) atop Mauna Kea. 

There have been several previous variability stud- 
ies of M 33. The first v ariab le stars in this galaxy were 
found bv lDuncanl Jl922f) and I Wolf I il923h . lHubtu7l lll926ri 
then found Cepheids in M33, thereby showing that this 
nebula is actually a gala xy at a distance of 263 kpc. 
iHubble fc Sandaee I dl953l) investigated a handful of the 
brightest variables in M33 and discovered the Hubble- 
Sandage class (n ow called Luminous Blue V ariables). Fol- 
lowing this work Ivan den Bergh et al.l lll975l) conducted a 
survey to find the brightest variables in M33 discovering 38 
new variables. The last major photographic variabil i ty sur - 
vey was conducted by iKinman. Mould and Wood I lll987t) . 
Most recently the DIRECT project conducted the first CCD 
variability survey of M33 discovering more than a thousand 
variables (Macri et al. 2001; Mochejska et al. 2001a, 2001b) 
in the centre of this galaxy. M33 has also been the target of 
several single (or few) epoch photometric surveys, with the 
survey by Massey et al. (2006, hereafter M06) as the most 
recent example. Our present survey has two improvements 
for finding variables over previous surveys of M33. The first 
is that by using a 3.6 metre telescope we go deeper than any 
previous survey. The second is that we are able to cover a 
much wider field (essentially the entire galaxy) than previ- 
ously possible for CCD surveys. As a result we have detected 
36709 variable sources which represents an increase by more 
than a factor of 20 in the number of known variables in the 
field of M33. 

In the following section we discuss the observations. In 
§3 we describe the data reduction and we present the cat- 
alogue of variable sources in §4. We describe a few of the 
ensemble properties of these variables in §5 and conclude 
with a few summary remarks in §6. 



2 OBSERVATIONS 

The data were obtained using the Queue Service Observing 
mode at the CFHT on 27 separate nights between August 
2003, and January 2005. We used the g' , r' and i' Sloan 
filters, obtaining 36 images centred on the galaxy with each 
filter. The seeing varied from 0.6 — 1" in £'. For the 2003 
observations we used an exposure time of 530 s for g , and 
660 s for r' and i' . In the 2004 season we switched to 480 
s for g' and 600 s for r' and i' . Table Q lists the dates of 
observations, the number of images obtained in each filter 
on these dates, and the median seeing for each filter. 

The MegaCam instrument is a wide-field mosaic imager 
consisting of 36 2048 x 4612 pixel CCDs. When used at the 
MegaPrime focus one obtains a 1 degree x 1 degree field of 
view with a sampling of 0'.' 187 per pixel. An example mosaic 
g' image is shown in fig.0 for comparison we also show the 
fields covered by DIRECT and the fields observed recently 
by M06. It is plain to see that with a single pointing it is 
possible to cover essentially the entire galaxy. 



Table 1. Observations of M33 



Date 


id 

9 


r 


i' FWHM g / 6 


"T~> T t 7TTTV If 

b WHM,./ 


"T~> T t 7TTTV If 

t WHM;/ 


2003-08-22 


1 


1 


1 


0.72 


0.70 


0.61 


2003-08-24 


2 


2 


2 


0.93 


0.91 


0.78 


2003-08-30 


1 


1 


1 


0.95 


0.82 


0.69 


2003-08-31 


2 


2 


2 


0.92 


0.90 


0.72 


2003-09-03 


1 


1 


1 


0.81 


0.75 


0.64 


2003-09-21 


1 


1 


1 


0.95 


0.84 


0.71 


2003-09-23 


1 


1 


1 


0.81 


0.79 


0.76 


2003-09-26 


1 


1 


1 


0.76 


0.67 


0.68 


2003-09-27 


1 


1 


1 


0.89 


0.85 


1.03 


2003-10-19 


4 


4 


4 


1.00 


0.92 


0.84 


2003-10-28 


1 


1 


1 


0.86 


0.85 


0.77 


2003-11-20 


1 


1 


1 


1.03 


1.14 


0.97 


2003-12-18 


1 


1 


1 


0.95 


1.06 


0.87 


2004-08-23 


1 


1 


1 


0.95 


0.99 


0.84 


2004-09-07 


1 


1 


1 


0.87 


0.85 


0.72 


2004-09-12 


1 


1 


1 


1.06 


0.79 


0.61 


2004-09-13 


1 


1 


1 


1.27 


1.20 


0.98 


2004-10-08 


1 


2 


3 


0.89 


0.99 


0.82 


2004-10-09 


1 


1 


1 


0.76 


0.75 


0.83 


2004-10-18 


1 


1 


1 


1.19 


1.20 


0.95 


2004-11-03 


1 


1 


1 


0.99 


0.82 


0.69 


2004-11-07 


3 


1 


3 


0.98 


1.00 


0.83 


2004-11-08 


1 


1 


1 


0.98 


0.89 


0.84 


2004-11-13 


1 


1 


1 


1.21 


0.96 


0.78 


2004-11-19 


1 


1 


1 


0.88 


0.85 


0.67 


2004-12-06 


1 


1 


1 


0.97 


0.83 


0.75 


2005-01-16 


1 


1 


1 


1.20 


1.17 


0.99 


a Number of 


images 


in 


this filter aquired 


on this date. 





Median FWHM in this filter 



3 DATA REDUCTION 

The preliminary mosaic CCD calibrations are performed as 
part of the CFHT Queue Service Observing mode. The im- 
ages are automatically processed through the Elixir pipeline 
which performs bias, dark, flat-field and fringe corrections, 
merges amplifiers, and provides an estimate of the photo- 
metric calibration to Sloan standards. 

The data reduction procedure consists of two distinct 
steps: in the first step we use image subtraction to iden- 
tify variable sources and obtain differential flux light curves, 
in the second step we perform point spread function (PSF) 
fitting to convert the differential flux light curves into mag- 
nitude light curves. The reductions were performed indepen- 
dently for each chip in each filter. 

To obtain light curves we use the image subtraction 
techniques due to Alard and Lupton (1998; Alard 2000) as 
implemented in the ISIS 2.1 package. This technique involves 
first registering the images to a master image. As a result 
we do not obtain light curves for any stars located on the 
chip gaps of the master image, we considered the simplicity 
of the standard ISIS scheme worth this sacrifice. We then 
combine a number of the best seeing images (we typically 
used 6 for each chip in each filter) into a master reference 
image. ISIS then solves for a transformation that matches 
the PSF, flux scale and background of the master reference 
and each image. The residuals between the transformed ref- 
erence and each image are then co-added (in absolute value) 
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Figure 1. A 36 CCD mosaic g' image of M33 obtained with the MegaPrime/MegaCam instrument on the 3.6 m. CFHT telescope. The 
field of view is f degree x 1 degree, north is up and east is to the left in this image. The smaller squares over the centre of the galaxy 
show the field of view of the DIRECT variability survey, while the larger squares show the coverage of the M06 photometric survey. The 
numbers for the chips are as used in the catalogue. 
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and searched for significant residual point sources which are 
identified as variables. A portion of the result from co-adding 
the residual images for chip 14 (as labelled in fig. in i is 
shown in fig. [5] together with the corresponding portion of 
the reference image. Each of the residual images has been 
divided by the square root of its original image before co- 
adding their absolute values to create this image (this corre- 
sponds to the "var" image produced by the findvar routine 
in ISIS). As a result, the value of each pixel is roughly pro- 
portional to its significance of variability. We also plot the 
variable sources from our catalogue in the image as open cir- 
cles with the radius of the circle proportional to log(Xjv_i) 
of the corresponding light curve where 

N „ 

i=l 

m,i is the magnitude at observation i, m is the average mag- 
nitude of the light curve, and <7; is the magnitude uncer- 
tainty for observation i. 

A number of the variable sources detected in this fash- 
ion are the results of artifacts in the data including sub- 
tle bad columns, slight jumps at amplifier boundaries, and 
diffraction spikes, bleed columns and reflections from par- 
ticularly bright stars. We attempt to eliminate the latter by 
growing masks over bright regions connected to saturated 
pixels. This is not 100% effective, and a few particularly 
pernicious regions are eliminated by hand. To eliminate ar- 
tifacts due to bad columns/rows we look for and remove 
spikes in the histograms of variables vs. column or row num- 
ber. In order to provide an unbiased estimate of the fraction 
of stars that are variable, we keep track of the eliminated 
regions and apply the same cuts to the point source lists 
determined below. We note that even with these cleaning 
procedures a number of artifacts may still slip into the cat- 
alogue. 

Once the variables have been identified we obtain their 
differential flux light curves by performing PSF fitting pho- 
tometry on the residual images. We obtained light curves for 
all variable sources that passed the above cleaning proce- 
dures as well as a number of randomly selected non-variable 
sources in each chip that can be used to measure the photo- 
metric uncertainty and determine the significance of detec- 
tion for the variables. 

In order to convert the differential flux light curves 
from image subtraction into magnitudes we perform PSF 
fitting photometry on the master reference images using 
the DAOPHOT/ALLSTAR package (Stetson 1987, 1991). 
Due to the large amount of data we ran the procedure in 
batch mode, relying on the the automated "pick" routine 
in DAOPHOT to choose PSF stars. We perform the pho- 
tometry in two passes, in the second pass we add in stars 
that are found on the image after subtracting off stars fit- 
ted in the first pass. For each chip we determine an optimal 
aperture correction radius to use in conjuction with ISIS 
that provides a minimum scatter among bright stars. We 
also determined an aperture correction at a fixed radius (13 
pixels, or roughly 3 times the FWHM) to allow for a mean- 
ingful comparison of magnitudes from separate chips. We 
performed the latter correction directly on the model PSF 
to avoid contamination due to crowding. Note that the flux 



scale for each master reference chip in a given filter comes 
from the same image. 

After performing PSF photometry we match the vari- 
able sources detected by ISIS to the point source catalogue 
from DAOPHOT. We only include variable sources that 
match within 2 pixels (0'.'37) to a point source. We use the 
DAOPHOT positions in measuring the light curves. As a 
final cleaning step we re-scale the magnitude uncertainties 
of all light curves on a given chip so that Xn-i as defined in 
eq.Qfor the constant light curves scatters about unity and 
does not correlate with magnitude. We then measure the 
standard deviation of Xn-i f° r the constant light curves 
and require that a variable source must have a light curve 
with Xn-i that is more than 5 standard deviations away 
from unity. Following the cleaning procedures we identified 
7640 variables sources in g , 15770 in r', and 27961 in i'. 

The final catalogue consists of 36709 distinct variable 
sources. We provide g' light curves for 26342 of these sources, 
r' light curves for 30386 and i! light curves for 34320. For 
a given variable we provide a light curve in every filter in 
which the object was detected by DAOPHOT. To match 
the point source catalogues from separate filters we used 
the ISIS interpolation routine to transform the g' and r' 
image coordinates to i' coordinates and then matched the 
catalogues with a 2 pixel (0'.'37) matching radius. There are 
201 close pairs (two objects separated by less than / .'3) in 
the final catalogue, we note in the comments section of the 
catalogue if a given source is a member of a pair. 

To provide approximate calibration to the Sloan stan- 
dard filters we used the photometric calibration provided by 
the Elixir routine. We will provide a better calibration with 
the full point source catalogue in a future publication. In im- 
plementing the Elixir calibration we only use the zero-point 
terms and ignore the colour terms. We do this for simplic- 
ity in applying the future calibration. Please note that the 
magnitudes provided in this catalogue should at present be 
regarded as instrumental magnitudes, they may differ from 
the standard system by as much as a few tenths of a magni- 
tude, depending upon the star's colour. We felt that the data 
are of sufficient value that it is worth releasing our variable 
stars catalogue, even though the photometric calibration is 
not definitive. 

We used the point source catalogue from M06 to deter- 
mine the astrometric solution for each chip. T he astrometry 
for M 06 is based on the USNO-B1.0 catalogue iMonet et alJ 
2003). For the few chips that do not overlap with M06 
we match directly to the USNO-B1.0 catalogue. The errors 
in this match were typically less than 0'.'2. For stars with 
R < 20 the median matching error to M06 is 0'.'13. 



4 CATALOGUE OF VARIABLE POINT 
SOURCES 

The catalogue of variable point sources and the g , r' and 
i! instrumental magnitude light curves are available on the 
web 1 , the full catalogue will also be included in the electronic 
version of the journal, in tables I2T11 we show a few rows to 



1 http: / /www. astro. livjm.ac.uk/~dfb/M33/ 
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Figure 2. Left: A l.'5xl.'5 portion of the i' reference image for Chip 14. Right: The corresponding "var" image from ISIS produced 
by co-adding in absolute value the residual images after rescaling by the square root of the original images. The circles show sources 
included in the catalogue of variables with the radius of the circle proportional to log(x^ r _ 1 ) of the light curve for that source. A few 
very bright regions on the "var" image and reference image are saturated stars that have been masked. 



Table 2. Catalogue of Variable Point Sources in M33 - 1. Chip numbers are as shown in fig. g' , r' and i' are the magnitudes 
measured on the reference image. The listed errors on these measurements are the formal uncertainties returned by DAOPHOT. 
g' flag, r' flag and i' flag denote whether or not the object was detected as a variable source in this filter. 



ID 


Chip 


«2000 Q 




glCCL 


g' err 


r> 


r" 'err 


i> 


i 1 err 


g'flag e 


r'flag 


i' 'flag 


10001 


1 


01:35:39.38 


31:00:51.2 


23.519 


0.023 


22.389 


0.009 


20.848 


0.004 





1 


1 


10002 


1 


01:35:39.79 


30:57:15.2 


24.185 


0.030 


22.093 


0.009 


21.137 


0.008 





1 


1 


10003 


1 


01:35:40.29 


31:05:36.1 


24.577 


0.026 


22.386 


0.006 


21.396 


0.005 





1 


1 


10004 


1 


01:35:43.17 


31:02:39.9 


22.924 


0.009 


99.999 


99.999 


20.369 


0.006 








1 


10005 


1 


01:35:44.04 


31:05:14.0 


99.999 


99.999 


23.443 


0.021 


21.158 


0.005 








1 


10006 


1 


01:35:44.95 


31:06:31.5 


23.878 


0.015 


22.522 


0.007 


20.769 


0.005 








1 


10007 


1 


01:35:47.46 


31:05:11.8 


99.999 


99.999 


99.999 


99.999 


23.394 


0.035 








1 


10008 


1 


01:35:47.57 


31:02:42.6 


24.646 


0.296 


99.999 


99.999 


99.999 


99.999 


1 








10009 


1 


01:35:47.64 


31:00:19.2 


23.927 


0.024 


22.568 


0.009 


20.820 


0.006 





1 


1 


10010 


1 


01:35:48.38 


31:06:15.0 


21.148 


0.005 


20.944 


0.005 


20.551 


0.004 


1 









a J2000, hours:minutes:seconds 

J2000, degrees:minutes:seconds 
c magnitudes 

d Value is 99.999 for no data. 

e Value is 1 if the object was detected as variable in this filter, otherwise. 



illustrate the available data. We split the catalogue here into 
four separate tables so that they each fit on a single page. 



4.1 Comparison with Other Catalogues 

The catalogue includes matches to both the variable star cat- 
alogues from DIRECT and the point source catalogue from 
M06. In fig. we compare the r' instrumental photometry 
from our constant stars with the R photometry from M06 
for chips 3 and 14. The transformation is strongly colour 
dependent, and varies from chip to chip. For each chip with 
at least 15 matches we fit a transformation of the form 

R = r' + aQ + al(g' -r') + a2(r -i'). (2) 



Table lists the resulting parameters of this fit along with 
the standard deviation after applying this fit to stars with 
R < 20. We take the standard deviation after applying an 
iterative 3er clipping to remove mismatches or stars that 
were resolved into different numbers of point sources in the 
two catalogues. We find that the scatter for bright stars 
is below 0.1 mag for all of the chips and below 0.03 mag 
for many of the chips that do not cover the centre of the 
galaxy. Though we do not apply these transformations to 
our catalogue, we provide them as a preliminary calibration. 

We also compared our catalogue with the lists of vari- 
ables found by the DIRECT project. We find matches within 
l'.'O to 1092 out of 1430 distinct variable sources. Of the 338 
variables that we do not match with, ~ 110 lie on or near 
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Table 3. Catalogue of Variable Point Sources in M33 - 2. The minimum and maximum magnitudes in each filter are the observed 
minimum and maximum magnitudes, no attempt has been made to cull out-lier points in determining these values, g'avg, r'avg 
and i'avg refer to the flux-averaged magnitude of the light curves, g'rms, r'rms and i'rms are the magnitude RMS values for 
the light curves, again determined without culling outliers. 



ID 


g'max ab 


g'min 


g'avg 


g'rms 


r'max 


r'min 


r'avg 


r'rms 


i'max 


i'min 


i'avg 


i'rms 


10001 


24.134 


23.131 


23.596 


0.246 


22.567 


21.907 


22.282 


0.185 


21.011 


20.510 


20.786 


0.132 


10002 


25.506 


23.962 


24.393 


0.366 


22.619 


22.018 


22.264 


0.198 


21.503 


21.020 


21.205 


0.146 


10003 


25.889 


24.371 


24.772 


0.387 


22.841 


22.295 


22.591 


0.172 


21.699 


21.339 


21.539 


0.117 


10004 


23.280 


22.764 


22.992 


0.141 


99.999 


99.999 


99.999 


99.999 


20.499 


20.325 


20.412 


0.051 


10005 


99.999 


99.999 


99.999 


99.999 


24.423 


22.976 


23.472 


0.274 


21.317 


20.851 


21.074 


0.115 


10006 


24.505 


23.446 


23.859 


0.253 


23.069 


22.091 


22.498 


0.232 


21.062 


20.569 


20.763 


0.135 


10007 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


24.384 


22.591 


23.142 


0.438 


10008 


27.074 


22.169 


24.352 


1.095 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


99.999 


10009 


24.796 


23.755 


24.098 


0.317 


23.480 


22.506 


22.845 


0.273 


21.429 


20.774 


20.996 


0.195 


10010 


21.331 


20.951 


21.192 


0.079 


21.077 


20.810 


20.984 


0.063 


20.681 


20.495 


20.591 


0.043 



a magnitudes 

b Value is 99.999 for no data. 

Table 4. Catalogue of Variable Point Sources in M33 - 3. IDm is the ID of the source in the M06 catalogue. Columns beginning 
with D33 are the IDs of the source in the DIRECT catalogues. Comments denote pairs of objects that are within 0''3 of one 
another, matches to Chandra sources (Grimm et al. 2005), or matches to H-cr emission line sources (Calzetti et al. 1995). Rows 
containing matches to M06 or DIRECT, or comments have been selected to illustrate the content of the catalogue. 



ID 


ID M a 




D33AID fe 


D33BID C D33ABID d 


Comments 


30073 


J013440.10- 


-305905.7 






Close Pair 30074 5 = 0.26; 


30074 










Close Pair 30073 6 = 0.26; 


30158 










Chandra J013444.6+305535 6 = 0.50; 


40584 


J013424.17- 


-305631.0 






Close Pair 40586 5 = 0.13; 


40586 










Close Pair 40584 5 = 0.13; 


130112 


J013407.23- 


h 304158.8 




D33J013407.2+304159.3 




130113 


J013407.21- 


^304635.0 


D33J013407.3+304635.5 


D33J013407.3+304635.5 




130349 






D33J013408.5+304430.6 


D33J013408.5+304430.6 




130408 


J013408.72- 


L 304543.1 


D33J013408.8+304543.5 


D33J013408.8+304543.5 




130504 






D33J013409.3+304238.6 


D33J013409.3+304238.6 





a Massey et al. (2006) 
6 Mochejska et al. (2001a) 
c Mochejska et al. (2001b) 
d Macri et al. (2001) 

a chip gap or near an image artifact. An inspection of the 
spatial distribution of DIRECT sources that did not match 
to our catalogue revealed that many of these sources are 
located near the densest part of the image where crowding 
is particularly significant and thus the coordinates are un- 
certain. W e find that when c omparing to the DIRECT AB 
catalogue (iMacri et al.ll200ll) the median matching error is 
O'.'ll, while for the DIRECT A and B catalogues (Mochejska 
et al. 2001a, 2001b) the matching error is 0'.'32 and 0'.'34 re- 
spectively. 

We have also matched our data to a catalo gue of Chan- 
dra X-Ray sources in M33 iGrimm et al]|2005 [). and a cata - 
logue of H-a line emission candidates~ ilcalzetti et al.lll995t) . 
We find that out of 261 distinct Chandra sources, 60 have 
an optical variable counterpart, where we used a matching 
radius of 1" plus the uncertainty in the Chandra position. 
For the match to Calzetti et al. (2005) we find that 58 out 



of the 153 sources have a counterpart in our catalogue when 
using a matching radius of Qf.'5. The IDs for these matches 
are provided in the comments column of the catalogue. 



5 RESULTS 

5.1 Colour-Magnitude Diagrams 

A colour-magnitude diagram (CMD) is perhaps the simplest 
plot one can make that reveals information about the classes 
of variables present in a catalogue. In fig0|we present g' — r' 
vs. r' and r' — i' vs. i' CMDs with and without the variables. 

In fig |S| we show the fraction of stars on the CMDs 
that were detected as variables. While biased, this plot al- 
lows us to identify several striking groups of variables. Two 
features that are immediately apparent are the Cepheid in- 
stability strip (IS) located between the main sequence and 
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Figure 3. A comparison between our photometry (g' ,r',i') and the photometry from M06 (Rmo6) f° r chips 3 (a-c) and 14 (d-f). There 
is a definite colour dependence in the transformation from r' to Rmo6 > plots (a) and (d) show the dependence of this transformation on 
g' —r' and plots (b) and (e) show the dependence on r' — i' . The bottom two plots (c,f) show the scatter after applying the transformation, 
here R' represents our transformed photometry. This transformation is preliminary, we will provide a definitive photometric calibration 
with the complete point source catalogue in a future publication. 



giant branches and the red long period variables (LPVs) 
at g' — r' > 1.25, r' > 20. Not unexpectedly there is 
a sizeable group of variables along the main sequence at 
—0.5 < r'—i' < 0, r' > 20, though the fraction of stars in this 
region that are variable does not increase dramatically. We 
also identify two groups of variables near the bright end of 
the main sequence and giant branches — at r' ~ 18 to 19, and 
g' — r' ~ -0.2 and +1.4, respectively. These we refer to tenta- 
tively as variable blue and red supergiant stars; they may be 
related to the a-Cygni type variables which are semi-regular 
pulsating supergiants (see for example Sterken 1996). The 



majority of variable candidates blueward of g' — r' = —0.5 
or r' — i' = —0.5 likely have incorrect magnitudes in at least 
one filter since they occur predominately in, and uniformly 
across, a handful of chips near the centre of the galaxy where 
magnitude measurements and matching between filters is 
the most uncertain. 

To find the blue and red edges of the IS, and the blue 
edge of the LPV region, we plot the fraction of stars that are 
detected as variable with RMS > 0.1 mag as a function of 
g' — r' and r' — i' colours in fig|S| We determine this fraction 
for stars with r',i' between 19.0 and 22.5 in bins of width 
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Figure 4. g' — r' vs. r' and r' — i' vs. %' CMDs. We have adopted the Hess technique of using intensity scaling proportional to the 
logarithm of the number of stars located at a given colour and magnitude. In the bottom two figures we use points to show the location 
of all variables in our catalogue on these CMDs. The preponderance of red variables along the giant branch and the location of the 
Cepheid instability strip between the main sequence and giant branch are apparent. 



0.5 mag. The location of the IS is visible in each of these 
distributions as a peak located between < g' — r' < 1 and 
< r' — i' < 0.5. Redward of the IS the fraction of stars 
that are variable rises until it reaches a break at around 
g — r' ~ 1.5 after which the rise is slower. We fit a simple 
model to the observed distribution of the form 



E 



(§)' 



+ Ix J 



+ 



(3) 



where x is exp(g'— r') or exp(r'— i'). The Ix J term is used to 
fit the distribution blueward of ~ 0, the first broken power- 
law is used to fit the IS and the second broken power-law 



is used to fit the LPV ridge. For the r' — i' distributions 
between 19.5 < i! < 21.5 we include a third broken power- 
law to account for the observed bump in the LPV region, 
particularly noticeable for the i' = 19.75 and i! = 20.25 
bins. Using this simple model we can define four interesting 
locations. The first is the peak in the IS given by (g' — 
r') pea k ~ hi(-B) (similarly for r' — i'). The second and third 
are the blue and red edges of the IS, which we will take 
to be the colours at which the fraction of stars that are 
variable falls to a tenth its value at the peak, and which 
we approximate by ( g P" me - r ') ISB = ln(B) + £ ln(20) and 
(g'-r')isit = ln(B) + i ln(20). The final location of interest 
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Table 5. Parameters for Photometric Transformation To R 
(M06). The parameters are denned in eq.|5] Only chips with more 
than 15 stars in common with M06 are listed. The scatter is the 
standard deviation after an iterative 3cr clipping for stars with 
R < 20. 



Chip ao ai a,2 Scatter 

0.041 
0.030 
0.042 
0.027 
0.085 
0.020 
0.093 
0.034 
0.045 
0.058 
0.034 
0.067 
0.044 
0.021 
0.024 
0.045 
0.065 
0.032 
0.024 
0.025 
0.026 
0.030 
0.019 
0.023 
0.018 
0.061 
0.027 



is the blue edge of the LPV region which we take to be the 
location of the first break redward of the IS and is given by 
(<?' - r') LPV = ln(F). 

In fig|7|we plot the measured positions of the peak in 
the IS, its blue and red edges, and the blue edge of the LPV 
region. We fit a line through the blue and red edges of the 
Cepheid IS and the blue edge of the LPV region and find 
that these locations are given approximately by the relations 



r'iSB = 


-10.64(±0.29)(ff' 


-r) +23.58(±0.86) 


(4) 


r'isR = 


-7.13(±0.15)(g' - 


r) + 26.27(±0.71) 


(5) 


r LPV — 


-18.32(±0.57)(fl' 


- r) + 45.93(±1.55) 


(6) 


i'lSB = 


-10.36(±0.30)(r' - 


- i') + 20.38(±0.84) 


(7) 


i'lSR = 


-19.94(±0.89)(r' - 


- i') + 27.72(±1.54) 


(8) 


i'lpv = 


-2.12(±0.01)(r' - 









+ 22.23(±0.11) 


i' < 20.75 


(9) 


■i 

% LPV — 


2.51(±0.02)(r' - i 


) 






+ 18.97(±0.17) 


i' > 20.75 


(10) 



Using these relations we show the approximate locations of 
IS and LPV stars in fig|7| We also identify regions around 
the red and blue supergiant branches that have an increased 
fraction of variability. 



5.2 Light Curves 

As a demonstration of our photometry we plot the g',r',i' 
light curves for several example variables from each of the 
four regions identified in figQ In fig|S]we plot phased differ- 
ential magnitude light curves for several Cepheids sorted by 
increasing period. We have set the faintest point in each light 
curve to have magnitude zero so that variations in colour as 
well as luminosity are apparent. In fig [5] we plot the light 
curves for several LPVs. We do not phase these light curves 
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Figure 6. The fraction of stars that are variable with RMS > 0.1 mag. as a function of colour (</ — r' left, r 1 — %' right) in magnitude 
bins of width 0.5 mag ranging from 19.25 to 22.25 (r' left, %' right). The error bars assume Poisson statistics. The solid lines show model 
fits to the observed distributions. The model is composed of a series of broken power-laws as described in §4.1. It is used to identify the 
peak of the Cepheid instability strip (between < g' — r' < 1 or < r' — i' < 0.5), its extent (which we take to be the colour at which 
the fraction of stars that are variable with RMS > 0.1 mag falls to a tenth of the peak value) and the blue edge of the long period 
variable population (which we take to be the location of the first power-law break redward of the Cepheid instability strip). 



or re-scale them. In fig llUl we plot differential magnitude light 
curves for several variable blue supergiant candidates, and 
in fie llll we plot differential magnitude light curves for sev- 
eral variable red supergiant candidates. The variable blue 
supergiants show variations on substantially shorter time- 
scales than the red variables, they also tend to have lower 
amplitudes than the red ones. 



5.3 Spatial Distributions 

In figs. 1121151 we show the spatial distributions of the vari- 
ables identified in fig. |7| As expected the long period vari- 
ables appear to trace the diffuse disk population, while the 
other classes of variables show a correlation with the spiral 
arms. Note that the variable red supergiant stars also ap- 
pear to be correlated with the spiral arms, this suggests that 
many of these variables are in fact red supergiant stars in 
M33 rather than foreground giants. 
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Figure 7. Same as fig. [S] Here we show classes of variables that we will use for further investigation. The solid lines bound the Cephcid 
region, the short-dashed lines bound the LPV region, and the long dashed lines bound the variable blue and red supergiants. We also 
plot the location of the Cepheid IS peak as found in fig. fusing triangles, its blue and red edges using squares, and the blue edge of the 
long period variable region using circles. 



5.4 Period-Magnitude Relation 

In fig. 1161 we show a preliminary period-magnitude dia- 
gram for the variables that lie in the Cepheid IS and have 
pulsation-like light curves. The magnitude plotted is the 
mean r' magnitude determined from fitting a Fourier se- 
ries to the light curves. We show this diagram to illustrate 
the quality and abundance of data available for Cepheids. 
We will conduct an analysis of the Cepheid population of 
M33 in a future contribution. 
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6 SUMMARY 

We have identified more than 36,000 variable sources in M33 
and are releasing photometry for all of them with multi-filter 
data for many of them. This makes M33 now comparable to 
either of the Magellanic cloud galaxies in terms of the sheer 
number of variable stars that are known. These data al- 
ready make possible a number of interesting projects. Stud- 
ies of the period-luminosity-colour relation for Cepheids in 
M33, and a test of its dependence on metallicity can be done 
with these data. This will be done when the photometry is 
properly calibrated. Investigations of the RMS or amplitude 
distributions of variables can be done using this relatively 
uniform data set. Studies of the eclipsing binary and ellip- 
soidal variable populations could be undertaken. And, with 
continued observations of this galaxy, it will be possible to 
study the LPV population for a galaxy beyond the Magel- 
lanic clouds. 
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Figure 10. Example differential magnitude light curves for 6 variable blue supergiants. We use open triangles for the g' observations, 
filled circles for r', and Xs for i' . The light curves have been shifted so that the faintest observation in each filter is at magnitude zero. 
This scale makes colour variations as well as luminosity variations apparent. 
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Figure 11. Same as fig. I1UI here we show variable red supergiants. 
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Figure 12. The location of Cepheid variable stars are plotted on top of a g' image of M33. The Cepheids appear to show some correlation 
with the spiral arms. 
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Figure 13. Same as fig. 1121 here we show LPVs. The LPVs do not appear to be correlated with the spiral arms. 
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Figure 14. Same as fig. 1121 here we show variable blue supergiant stars which do appear to be correlated with the spiral arms. 
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Figure 15. Same as fig, 1121 here we show variable red supergiant stars which do appear to be correlated with the spiral arms. 
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Period [days] 

Figure 16. The period-magnitude relation for ~ 2500 variables located in the Cepheid IS that have pulsation-like light curves. The 
magnitude plotted is the mean r' magnitude of the variable determined by fitting a fourier series to the light curve. No effort has been 
made to correct for extinction. Cepheids pulsating in the fundamental and first overtone modes arc clearly visible on this diagram. 



